Hierarchical freezing in a lattice model 
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A certain two-dimensional lattice model with nearest and next-nearest neighbor interactions is 
known to have a limit-periodic ground state. We show that during a slow quench from the high 
temperature, disordered phase, the ground state emerges through an infinite sequence of phase 
transitions. We define appropriate order parameters and show that the transitions are related 
by renormalizations of the temperature scale. As the temperature is decreased, sublattices with 
increasingly large lattice constants become ordered. A rapid quench results in glass-like state due 
to kinetic barriers created by simultaneous freezing on sublattices with difi'erent lattice constants. 



A recent result in tiling theory [1| presents a new op- 
portunity for studying the development of long-range or- 
der in a system with a non-periodic ground state. A sin- 
gle hexagonal prototile has been shown to force a limit- 
periodic pattern (meaning a state made up of the union 
of an infinite number of periodic structures with lattice 
constants of ever-increasing sizes 0-0] ) ■ The existence of 
a local Hamiltonian with a limit-periodic ground state 
composed of a single unit repeated in different orien- 
tations suggests that it may be possible for solid state 
materials, colloidal systems, or perhaps even collections 
of macroscopic units to realize structures of this type. 
A similar type of complexity occurs in quasicrystals and 
has implications for electronic, photonic, elastic, and fric- 
tional properties. |5j Even in the absence of any physical 
examples of limit-periodic phases, however, the statis- 
tical mechanical properties of the new tiling model are 
remarkable and may serve as a conceptually useful inter- 
polation between crystalline and glassy behavior. 

In this paper, we study the spontaneous formation of 
the limit-periodic structure. Working with a substan- 
tially more complex tiling model with square symmetry, 
Mi§kisz showed that a series of partially ordered equi- 
librium states exist at nonzero temperatures in a system 
with a limit-periodic ground state Here, we explic- 
itly define order parameters for the transitions, present 
numerical evidence for ordering at slow cooling rates 
through an infinite sequence of phase transitions, and 
present a scaling theory of the transition hierarchy. We 
also show that rapid quenching produces a state with 
frozen disorder resulting from the competition between 
two or more levels of the hierarchy. 

Our model Hamiltonian is based on rules for placing 
hexagonal tiles of the type shown in Fig. [Ha) on a close- 
packed lattice. Each tile may be placed in any of the 
twelve orientations obtainable by rotations by 7r/3 and re- 
flection. For each nearest neighbor pair, we assign energy 
zero if the black stripe is continuous across their shared 
boundary and ei otherwise. For each next-nearest neigh- 
bor pair, the energy is zero if the flags in their closest cor- 
ners point in the same direction and €2 otherwise. Ref. [l] 
shows that for positive ei and 62, the ground state of an 
infinite system is a zero-energy structure that contains 



black triangles with arbitrarily large side lengths. We 
emphasize that the system is completely homogeneous, 
being a lattice of identical units. 



STAGGERED TETRAHEDRAL ORDER 

The key to explaining the behavior of the lattice model 
is to focus on the patterns of (truncated) black triangles 
of different sizes in Fig. [Ijc). We refer to the smallest 
triangles as making up level 1, the next smallest as level 
2, etc., and note that the edge of a triangle at level n is 
a straight black stripe crossing /c„ — 1 tiles, where fc„ = 
2"~^. At each level we have a periodic arrangement of 
triangles centered on the vertices of a honeycomb lattice. 

The level- 1 pattern may form in four different ways. 
One example is shown in Fig. [21 in which the tiles marked 
A do not contribute at all to the triangles. The other 
three are translations of this one, each with a different 
non-contributing sublattice (B, C, or D). We associate 
with each tile j a "staggered tetrahedral spin" vector 
(Ti j- = ex J where X indicates one of the four vertices of 
a reference tetrahedron. (See Fig. [2j) The spin of a tile is 
determined both by the orientation of the diameter join- 
ing its two black triangle corners and by the sublattice 
that it belongs to, according to the map shown at right 
in Fig. [2] For example, a tile with corners aligned verti- 
cally and sitting on the D sublattice is assigned cti = ec- 
Note that specifying <Ti j does not completely specify the 
orientation of tile j. There are four consistent choices, 
corresponding to the two possible locations of the long 
black stripe and two possible orientations of the flag bar 
perpendicular to that stripe. Note also that for any given 
tile, <Ti can take only three of the four possible values. 

We define the average total spin cri tot = Tlj ""ij j 
where N is the number of tiles in the system. The system 
exhibits tetrahedral symmetry in the following sense: for 
each configuration with a given CTi^tot, there is another 
with identical energy having ct'i ^ot related to <Ti_tot by an 
operation in the 24-element tetrahedral group T^. The 
mapping from operations on the lattice to elements of Td 
is given in Table [Ij 

Consider now the level-2 triangles. If the level-1 struc- 
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Lattice operation Td operation on (f) 

Rotation by 2n/3 Rotation by 2n/3 

about center of X about ex 



Reflection through edge Reflection through 

shared by X and Y (ex,ey) plane 



Translations taking — 5- Rotations by vr 
X sublattice to Y about ex + ey 



Rotation by 2n/3 — >■ Rotary inversion 
followed by reflection 



6T = 0.01 after every Nt steps, where t controls the 
quench rate. 

For the rapid quench (see Fig. [3]) , <j)i reaches its sat- 
uration value but (p2 does not. The ordering on level 
2 ceases to increase once the level-3 transition temper- 
ature is reached, at which point interactions associated 
with level 3 prevent further ordering of level 2. For the 
slow quench, on the other hand. Fig. |4] shows a sequence 
of clear transitions. The number of transitions that can 
be observed is limited by the system size. At level 5 (not 
shown), the edge length of a black triangle is 15 tiles, so 
our 64 X 64 system is only 4 x 4 at level 5. 

Because the symmetry of the order parameter is the 
same as that of a 4-state Potts model, we expect the 
transition to fall in the same universality class, which 
has an order parameter exponent /3 = 1/12. Q Prelimi- 
nary numerical investigations are consistent with a very 
sharp second order transition, having not revealed any 
hysteresis in a cooling and heating cycle. None of the 
analysis below depends on the critical behavior or even 
on the order of the transition. 



TABLE I. Symmetry operations for the total staggered tetra- 
hedral spin. The left column specifles an operation on the 
2D tiling pattern, where X,Y £ {A, B ,C, D} each represent 
a tile in the corresponding sublattice of Fig. [51 The right col- 
umn specifies 3D operations on the order parameter in terms 
of the tetrahedral star of vectors ex, where X is the label 
shown on Fig. [5] 



ture is perfectly ordered, the corners of the level-2 trian- 
gles must come from the non-contributing sublattice at 
level 1. That sublattice (e.g., all of the A's in Fig. ^ 
is a precise copy of the original lattice (with distances 
scaled by a factor of 2), so we can define a new stag- 
gered tetrahedral spin cr2^tot = wjj^j where the 
sum runs only over the sublattice of interest. The con- 
struction can be repeated ad infinitum, with (t„ tot being 
well-defined if and only if the ordering on level n — 1 is 
sufficiently strong to unambiguously specify which sub- 
lattice will order at level n. 



FRUSTRATION AND ORDER IN QUENCHES 

We define an order parameter for each level of struc- 
ture: (pn = max[ex • f n.tot], where X € {A, B, C, D} and 
ex is a unit vector in the direction labeled "X" in Fig. [21 
We take the magnitude a of each spin to be 3/2 so that 
the maximum 0„ is unity for all n. The maximum oc- 
curs if and only if the lattice of small triangles is perfectly 
ordered, regardless of the orientations of the tiles on the 
X sublattice. Figs. [31 and [H show the behavior of 0„ with 
n — 1,2,3,4 for a rapid quench and for a slow quench. 
The simulations were done on rhombic domain with pe- 
riodic boundary conditions. Because the ground state is 
nonperiodic, there must be some mismatches at all times. 
For domains of size 2" x 2", the smallest possible number 
of mismatches is 4. The simulations employed a standard 
Metropolis algorithm in which a random tile is chosen at 
each step and one of the twelve orientations of that tile 
is chosen randomly as a possible move. T is lowered by 



SCALING THEORY 

The system is considered fully ordered at level n when 
all of the level-n triangle corners occur in the pattern 
of Fig. mjc). In the ordered pattern, at each level, fiags 
in the tiles surrounding a non-contributing tile also form 
triangular structures of the type shown at left in Fig. [2l 
The long black stripes and the flag bars perpendicular 
to them, however, need not be in their ground state ori- 
entations. Their role is only to mediate the interactions 
between the triangle corners. We will see below that if 
the triangle corners at all levels m < n are assumed to be 
perfectly ordered and immovable, the level-n dynamics 
are identical to the level- 1 dynamics at a rescaled tem- 
perature (and possibly a rescaled £2). We discuss first 
the special case ei — 62 = e. 

The partition function for the level-n system can be ex- 
pressed as follows. A configuration of the system can be 
specified by giving, for every edge of each hexagon, the 
location of the black stripe meeting that edge and, for 
every vertex, the orientation of the fiag at that vertex. A 
configuration is allowed if and only if the specification for 
every hexagon corresponds to some orientation of the tile 
decoration of Fig. [Ija) , independent of whether neighbor- 
ing tiles match properly. At level n, the contribution to 
the partition function from a given configuration of the 
triangle corners is a product of contributions from all of 
the nearest and next-nearest neighbor bonds between the 
corner tiles. Each of these bonds consists of fc„ — 1 tiles, 
with black stripes joining nearest neighbors of the level- 
n sublattice and flag bars joining next-nearest neighbors. 
Fig. [21 shows an example of each type of bond for the 
n = 2 case. The dashed grey lines show the two possible 
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locations of the stripe on a triangle edge joining the cor- 
ners that would be present on the A sublattice. The grey 
diagonal shows the two possibe orientations of the flag 
bar that forms a triangle edge for next-nearest neighbor 
interactions between two A tiles. 

The full partition function for the level-n system is 



ZniT) = n Cn.fc 



(1) 



configurations \bonds 



where b € {odd, even} represents the state of a given 
bond (mismatched or matched) in the given configura- 
tion. Each edge is effectively a ID Ising system with A;„ 
possible mismatches (see Fig. [5l^a)), which gives 



dd 



(l + e-^/^)'" + (l-e-/^)' 



(2) 



Now because the configuration sums are identical for 
all levels, the behavior of the system at level n at some 
temperature T„ will be identical to that for level 1 at tem- 
perature Ti if Cn,odd(T'„) = aCi,odd(T'i) andCn,ovon(r«) = 
« Ci.cvcn(T'i) for some normalization constant a. Elimi- 
nating a from these equations yields, after some algebra, 



tanh ( — |- I — 
or, equivalently, for all n. 



tanh 



2T„ 
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2T„ 
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From Eq. ([3]), the transition temperature for large n is 
1 



n log 2 — log 
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(5) 



One might worry that Eq. ([5]) will break down because 
may not be fully saturated at Tc-n+i- Nevertheless, as 
shown in Fig. [SI we obtain an excellent data collapse for 
several levels by plotting 0„(T„) as a function of Ti(T„). 
(Note: The x 's on the high temperature side of the tran- 
sition represent finite-size fiuctuations and correspond to 
projections onto different tetrahedral vectors.) The col- 
lapse shows that the residual disorder after the transi- 
tion at level n has a very small effect on the transition 
at level n -I- 1. For n = 1, the data for Fig. |4] show 
clearly that is very close to unity (~ 0.999) at the 
critical temperature for (f)2. The scaling argument takes 
this value to be exactly unity, in which case Eq. (|4]) leads 
to (pn(Tc-n+i) = <i>\iTc\'i) for all n. If the small deviation 
from unity causes the scaling to break down at large n, 
it appears that the difficulty will only emerge at lattice 
sizes too large to be probed computationally. 



The analysis is more complicated for the generic case 
£2 < ei, but the essential phenomenology appears to be 
the same. Let T'c;„(ei,e2) be the level-n transition tem- 
perature. Assuming that deviations discussed in the pre- 
vious paragraph can still be neglected, T^-^n must be an 
increasing function of each of its arguments. Thus we 
must havee2/ei < T'c;„(ei, e2)/Tc;„(ei, ei) < 1, where the 
first inequality follows from the fact that a simple rescal- 
ing of aU energies gives Tc;„(e2, £2) = (e2/ei)T'c;„(ei, ei). 



DISCUSSION 

The ground state of the tiling model is discussed in 
detail in Ref. [l| . The complexity of the structure is best 
revealed by the pattern of tile parities (the grey and white 
in Fig. [ijc)). We have shown here that the model also 
has remarkable properties at finite temperatures. First, 
it exhibits a hierarchy of distinct thermodynamic phases, 
with each successive one corresponding to formation of 
an ordered lattice with a lattice constant twice as large 
as the last. Second, each transition is unusual in that 
the low temperature phase leaves one quarter of the tiles 
as "rattlers" with undetermined orientations. Finally, 
the kinetics of ordering become frustrated if the quench 
rate is fast enough that the temperature drops below Tc;„ 
before c^n-i has reached a sufficiently high value. The 
latter effect is reminiscent of glass formation, where lower 
energy states can be reached through slower quenching 
but any nonzero quench rate eventually leads to trapping 
in a nonequilibrium configuration. Q The present model 
shows that this can happen in the context of a series of 
transitions that each establish true long-range order. 

Many questions remain concerning the precise nature 
of the staggered tetrahedral phase transition, the pos- 
sibility of reaching the ground state via growth from a 
small seed, the geometry and energetics of topological 
defects in this system, and the scaling of the time re- 
quired for ordering at each level. 

To date, we know of no examples of spontaneous emer- 
gence of limit-periodic order in a real system. The 
present work suggests that such systems may exist and 
exhibit novel properties. The hexagonal tiling consists of 
a single unit that may be realized as a cluster of atoms, 
a surface pattern on a colloidal particle, or even a micro- 
machined brick. It is also possible that such systems have 
already been observed but not recognized because they 
present as "poor" crystals in diffraction studies. (See 
Ref. [9] for a discussion of diffraction from a limit-periodic 
structure.) Systems of interest may be two-dimensional 
monolayers on a fiat substrate, surface reconstructions on 
crystals, or bulk phases. See Refs. [ij for a 3D tile that 
enforces the present structure through its shape alone by 
allowing direct contact between next-nearest neighbors. 

Mi§kisz's study based on an extension of Robin- 



son's set of six Wang tiles [11| indicates that a limit 
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periodic ground state with square symmetry should show 
a similar sequence of transitions. See [1^ and 13 1 for 
2D and 3D examples of limit-periodic tilings with square 
symmetry composed of only two types of tiles, which 
may be more amenable to the type of analysis performed 
above. 
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FIG. 1. The hexagonal prototile and its mirror image, (a) The two tiles are related by reflection about a vertical line, (b) For 
zero energy, adjacent tiles must form continuous black stripes and flag decorations at opposite ends of a tile edge (as indicated 
by the arrows) must point in the same direction, (c) A portion of an infinite tiling that has zero energy. 




FIG. 2. The sublattices employed in the definition of the order parameter and the mapping from tile orientations to spin 
vectors. The tiles of the A sublattice are labeled along with one tile each of the B, C, and D sublattices. The figure shows the 
pattern of level- 1 triangles formed when the non-contributing tiles are those of the A sublattice. For explanation of the cluster 
shown at left, the dashed lines, and grey bar, see text. 



FIG. 3. Right: The behavior of the order parameters for a rapid quench. Simulations were performed on a 64 x 64 rhombic 
domain for ei = £2 = 1 and t = 120. Left: The result of a rapid quench on an 8 x 8 lattice showing the high density of defects 
that persists at long times. There are no defects, however, in the level-1 structure, consistent with the high value of 01. 
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FIG. 4. Right: The behavior of the order parameters for a slow quench. Simulations were performed on a 64 x 64 rhombic 
domain for ei = £2 = 1 and r = 12 x 10^. Left: The result of a slow quench on an 8 x 8 lattice showing the minimum possible 
number of defects. 
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FIG. 5. (a) Matched and mismatched corner configurations, shown here for level-3 triangle edges, (b) Scaling collapse of Fig.U 
data. The deviation of the level-4 points on the right is due to the finite size of the system. 



